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Nomenclature 


A  =  nozzle  section 

h  =  Plank  constant 

Cp  =  constant  pressure  heat  capacity 

Cv  =  constant  volume  heat  capacity 

c  =  speed  of  light 

Ea  =  reaction  activation  energy 

e  =  electron  charge 

Kd  =  rate  coefficient  of  the  direct  reaction 

Keq  =  equilibrium  constant 

Kr  =  rate  coefficient  of  the  reverse  reaction 

Ma  =  Mach  number 

P0  =  reservoir  pressure 

R  =  universal  gas  constant 

Tn  =  reservoir  temperature 

Tv  =  vibrational  temperature 

v  =  vibrational  quantum  number 

a  =  Cv/R 

a/,  =  ammonia  dissociation  degree 

sY  =  vibrational  level  energy 

Y  =  Cp/Cv 

we  =  first  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

0):xe  =  second  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

(0:.ye  =  third  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

(0,,ze  =  fourth  order  spectroscopic  constant  to  calculate  vibrational  level  energies 


I.  Introduction 


ecently,  solar  thermal  propulsion  (STP)  engines  have  been  investigated  as  a  promising  technology  to  enhance 
microsatellite  maneuverability  and  enable  new  missions.  Solar  thermal  propulsion  takes  advantage  of 
concentrated  sunlight — typically  at  concentration  ratios  of  up  to  10000:1 — in  order  to  raise  a  heat  exchanger 
to  temperatures  of  several  thousand  degrees  K.  At  the  Surrey  Space  Centre1'2  (SCC),  a  small  scale  STP  engine  has 
been  designed  to  use  ammonia  as  propellant.  The  main  advantage  of  ammonia  is  that  it  is  a  storable  liquid  with  low 
molar  mass  and  easy  handling  characteristics  and  does  not  exhibit  the  serious  toxicity  and  flammability  risks  posed 
by  hydrazine. 

To  achieve  improved  performances,  the  ammonia  gas  is  heated  to  2000-2500  K,  which  should  result  in 
dissociation  to  N2  and  H2  molecules,  doubling  the  particle  density  and  decreasing  the  average  molecular  weight  of 
the  mixture.  This  improves  rocket  performance  by  raising  specific  impulse  (Isp),  which  varies  as  the  inverse  square 
root  of  molecular  weight.  The  problem  is  that  the  dissociation  kinetics  at  temperatures  of  interest  in  the  engine  are 
very  slow  and  therefore  the  complete  dissociation  of  ammonia  molecules  may  not  occur. 

In  this  paper  we  have  developed  a  kinetic  model  for  ammonia  mixtures  applied  to  the  nozzle  investigated  at 
SSC1'2.  It  is  worth  noting  that  we  have  introduced  state-to-state  kinetics  to  verify  the  presence  of  non-equilibrium 
vibrational  distributions  during  the  nozzle  expansion.  We  have  also  investigated  the  behaviors  of  the  engine  for 
higher  reservoir  temperatures  (3000K-5000K)  where  the  dissociation  of  ammonia  is  complete. 


II.  Model  Description 


A.  Nozzle  Geometry 

The  nozzle  used  in  the  thruster  has  a  conic  geometry  of  the  following  dimensions  (in  mm): 


throat  radius  0.35 

exit  radius  2.0 

exit-throat  distance  6.2 

inlet  radius  1.0 

inlet-throat  distance  1 .2 


These  data  have  been  provided  by  SSC  researchers.  The  converging  nozzle  is  quite  short  so  we  choose  to  include 
part  of  the  reservoir  in  the  calculation  through  a  connecting  cone  of  the  following  dimensions  (in  mm) 


reservoir  radius  1 1 .25 


connecting  cone  length  2.0 
The  nozzle  is  shown  in  fig  1 . 

The  fluid  dynamic  model  of  the  nozzle  flow  is  based  on  the  almost  one-dimensional  Euler  equations  for  a 
variable  section  duct.  The  Euler  equations  are  solved,  for  the  stationary  state,  together  with  the  perfect  gas  equation 
of  state  and  the  mass  continuity  for  each  species  including  the  reaction  contributions.  The  numerical  scheme  adopted 
is  described  in  ref. '. 

B.  Chemical  Model 

To  build  the  chemical  kinetic  models,  we  need  to  select  both  the  species  that  must  be  taken  into  account  and  their 
mutual  reactions.  In  the  initial  demonstration  laboratory  experiments,  the  pressure  ranges  from  1  to  4  atm  and  the 
temperature  is  less  than  1500  K.  For  operational  versions,  the  temperature  could  reach  the  value  of  2500  K.  The 
pressure  values  are  measured  directly,  but  the  gas  temperature  is  estimated  from  some  measurements  performed 
outside  the  reservoir.  This  results  in  a  large  uncertainty  on  the  actual  conditions  operating  inside  the  reservoir.  For 
these  temperatures  the  ammonia  decomposition  kinetics  are  quite  slow.  Surface  dissociation,  neglected  in  this 
approach,  should,  probably,  be  inserted  in  the  model. 

1.  Species 

We  assume  that  the  ammonia  mixture  consists  of  atoms  and  molecules  formed  in  the  NH3  decomposition.  We 
ignore  all  the  ions  because  the  temperatures  involved  are  quite  low.  We  introduce  the  vibrational  excitation  only  for 
N2  and  Eli  molecules,  which  have  stable  vibrational  excited  levels.  For  all  the  other  molecules,  we  disregard  the 
internal  levels  as  their  number  is  so  high  that  the  computational  load  becomes  unaffordable.  Moreover,  there  is  a 
lack  of  data  for  implementing  state-to-state  kinetics  for  NF1X  molecules.  A  multitemperature  model  could  be  a  valid 
approach  assuring  an  acceptable  computational  cost  nevertheless,  in  this  work,  we  consider  these  internal  states  in 
equilibrium  with  the  translational  degree  of  freedom  or  as  following  a  frozen  kinetics,  simulating  non-equilibrium 
expansion.  According  to  the  internal  state  model  selected  we  assign  the  Cp  and  y  (y  =  Cp/Cv)  coefficients  to  each 
molecule.  For  each  vibrational  mode  or  rotational  axis,  in  equilibrium  with  the  translational  degree  of  freedom,  we 
add  the  quantity  aR  into  the  expressions  of  Cp  and  Cv,  where  R  represents  the  universal  gas  constant  and  the  a 
parameter  is  set  to  1  for  a  vibrational  mode  and  to  Vi  for  a  rotational  axis.  This  approximation  is  valid  for  the  rigid 
rotor  and  for  the  infinite  harmonic  oscillator,  since  in  general  the  Cp  coefficient  depends  on  the  temperature.  For 


atomic  species  no  internal  degree  of  freedom  is  considered.  The  contribution  of  translational  degrees  to  the  molar 
specific  heat  is  2.5R.  Table  1  reports  the  species  accounted  and  their  relative  properties. 

For  H2  and  N,  we  estimate  the  energy  of  each  vibrational  level  (sv)  through  the  semi-empirical  formula4, 
function  of  the  vibrational  quantum  number  and  of  some  spectroscopic  constants  expressed  in  cm"1: 
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The  H2  spectroscopic  constants  have  been  obtained  through  the  fits  of  the  vibrational  level  energies  extracted  from 
the  BKMP  potential  energy  surface5  while  for  N2  we  used  data  from  ref.  6.  Table  2  summarizes  those  values  for  the 
two  species. 

2.  Processes 

We  implement  in  our  model  macroscopic  kinetics  for  NH3  depletion  based  on  the  following  processes: 


nh3+x  <-»  nh2+h+x 

(pl) 

NH3  +  X  NH  +  H2  +  X 

(p2) 

NH3  +  H  <-»  NH2  +  H2 

(p3) 

NH  +  H  +  X  NH2  +  X 

(p4) 

H2  +  N  +  X  <->  NH2  +  X 

(p5) 

NH2  +  H  <-4  NH  +  H, 

(p6) 

NH2  +  NH,  <-»  NH3  +  NH 

(p7) 

N  +  H  +  X-c^NH  +  X 

(p8) 

NH  +  H  H,  +  N 

(p9) 

NH  +  N  <->  N,  +  H 

(plO) 

NH  +  NH  H,  +  N, 

(pH) 

2H  +  N,  H,  +  N, 

(pl2) 

2H  +  H,  <->  H,  +  H, 

(p  13) 

2N  +  N,  -C4  N,  +  N, 

(p  14) 

2N  +  N  <->  N2  +  N 


(Pi  5) 

NH,  +  NH,  ->  N2H3  +  H2  (p  16) 

N2H3  +  X^NH2  +  NH  +  X  (pi  7) 

N2H3  +  H  — »  2NH2  (pi 8) 

N2H3  +  H  -»  NH3  +  NH  (pi 9) 


where  X  is  a  generic  component.  The  rate  coefficients  for  all  these  processes  are  expressed  in  the  Arrhenius  form4 


K  =  K0T  e 


6P-Ea/T 


(2) 


For  the  direct  reaction,  the  parameters  K0,  5  and  Ea  are  reported  in  Tab.  3.  The  rate  coefficients  of  these  processes 
have  been  taken  from  ref.  6-11.  In  this  model,  we  neglect  more  complex  molecules,  such  as  N2H4,  that  do  not  seem 
to  be  significant. 

Rates  for  the  reverse  processes  have  been  obtained  through  the  detailed  balance  principle4 


K„  = 


K, 


K, 


eq 


(3) 


where  K<,q  is  the  equilibrium  constant.  Equilibrium  constants  have  been  calculated,  from  partition  functions, 
following  the  classical  statistical  thermodynamic  theory12.  Subsequently  they  have  been  fitted  by  the  equation 


ln^eq)=K00+Kp 


KX)0 
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(4) 


We  consider  only  the  equilibrium  constants  for  the  generation  of  a  new  species  from  standard  ones  like  N2  and  Fl2  in 
nitrogen-hydrogen  compounds  formation.  Therefore  all  the  formation  reactions  can  be  synthesized  as 


n 

2 


N2  + 


m 

~2 


H2  <-»  NnHm 


(5) 


The  parameter  values  utilized  to  calculate  Kcq,  for  the  species  listed  in  table  1  are  reported  in  tab.  4.  From  these 
formation  equilibrium  constants  we  can  calculate  the  equilibrium  constants  relative  to  all  the  reactions  involving  the 
species  in  the  table. 

The  kinetic  model  has  been  improved  including  state-to-state  kinetics  for  the  diatomic  species  N2  and  H2, 
substituting  processes  P12-P15  with  the  state-selective  dissociation  and  adding  vibrational  relaxation  processes13"15 


N2(v)  +  N2(w)  <-»  N2(v-1)  +  N2(w+1) 

(p20) 

N2(v)  +  N2(w)  <-»  N2(v-1)  +  N2(w) 

(p21) 

N2(v)  +  N  <->  N2(v-Av)  +  N 

(p22) 

N2(v)  +  N2(w)  <-4  N2(v-1)  +  2N 

(p23) 

N2(v)  +  N2(w)  <->  2N  +  N2(w) 

(p24) 

N2(v)  +  N  2N  +  N 

(p25) 

H2(v)  +  Pl2(w)  <->  F12(v-1)  +  H2(w+1) 

(p26) 

H2(v)  +  Pl2(w)  <->  F12(v-1)  +  H2(w) 

(p27) 

H2(v)  +  Ho  H2(v-Av)  +  H 

(p28) 

H2(v)  +  H2(w)  H2(v-1)  +  2H 

(p29) 

H2(v)  +  H2(w)  2H  +  H2(w) 

(P30) 

H2(v)  +  H  <-*  2H  +  H 

(p31) 

All  the  other  processes  involve  only  the  ground  state  of  such  molecules. 

relaxation  of  nitrogen  molecules  due  to  hydrogen6'9: 

We  have  also  included  vibrational 

N2(v)  +  Pl2(w)  <-»  N2(v-1)  +  H2(w) 

(P32) 

N2(v)  +  H  <-4  N2(v-Av)  +  H 

(p33) 

N2(v)  +  Pl2(w)  <->  N2(v-1)  +  H2(w+1) 

(p34) 

III.  Results 


C.  Low  reservoir  temperature  (T<2500  K) 


In  this  paragraph  we  want  to  investigate  the  kinetic  behavior  of  the  ammonia  mixture  in  cases  close  to  the 
working  conditions  of  the  solar  thruster  studied  at  the  Surrey  Space  Centre1'2.  The  pressure  range  is  between  1  and  4 
bars  and  the  temperature  is  less  than  1500  K.  The  heater  temperature  is  not  well  known  because  it  has  been 
measured1,2  outside  the  chamber.  We  have  extended  our  calculations  to  temperature  values  below  2500  K,  which 
may  be  obtained  in  orbital  conditions. 

We  investigate  four  different  pressure  values  in  the  range  1-4  bars.  In  figure  2  we  report  the  concentration 
profiles  along  the  nozzle  obtained  for  reservoir  conditions  of  P=4  bars  and  T=1500  K,  and  introducing  the 
vibrational  kinetics  of  N2  and  H2.  It  is  interesting  to  observe  that  the  flow  is  practically  frozen  along  the  whole 
nozzle  except  for  some  minor  species  such  as  H  and  NH2  that  completely  disappear  at  the  nozzle  exit.  The  atomic 
nitrogen  molar  fraction  is  lower  than  10"10.  From  fig.  1  it  can  be  concluded  that  the  contribution  of  NHx  molecules 
to  the  flow  characteristics  is  negligible  because  the  concentration  is  very  small. 

Figures  3  and  4  describe  the  vibrational  distributions  of  hydrogen  and  nitrogen  for  the  conditions  of  figure  2.  We 
can  observe  that  the  nitrogen  distributions  are  frozen.  On  the  contrary,  while  H2  distributions,  for  v<3,  are  almost 
frozen,  their  tails  are  strongly  depleted,  resulting  in  non-Boltzmann  distributions.  The  different  behavior  of  the  two 
species  is  due  to  the  different  percentage  of  nitrogen  and  hydrogen  atoms;  hydrogen  atom  concentration  is  high 
enough  so  that  the  VT  collisions  involving  atoms  are  frequent  enough  to  cool  the  distribution  tails. 

The  distribution  behavior  brings  two  conclusions:  first  of  all  the  contribution  of  vibrational  kinetics  is  practically 
negligible  and  the  vibrational  temperature  can  be  considered  frozen,  at  least  for  diatomic  molecules.  As  a 
consequence,  in  the  following  results  we  have  used  only  the  macroscopic  model  freezing  the  vibrational  temperature 
of  all  the  molecules,  in  the  sense  that  the  contribution  of  the  vibrational  degrees  of  freedom  is  null.  This  results  in 
neglecting  the  vibrational  contribution  (“Vib”  in  Tab.  1)  in  the  calculation  of  Cp. 

We  observe  that  the  Ma  is  practically  independent  of  the  initial  temperature  and  pressure,  obtaining  at  the  nozzle 
exit  a  value  close  to  Ma=5A.  We  report  the  pressure  along  the  nozzle  in  figure  5.  It  is  quite  interesting  to  note  that 
the  exit  pressure  is  proportional  to  the  inlet  pressure.  These  two  results  occur  as  the  chemistry  weakly  affects  the 
macroscopic  quantities  during  the  expansion. 

Table  5  reports  the  mass  flow  rate  and  the  thrust  for  different  inlet  pressures  at  To=1500  K  or  To=2500  K.  Table 
6  on  the  other  hand  reports  experimental  results  of  mass  flow  rate  and  thrust  measured  on  the  MKI  solar  receiver  by 
Kennedy  and  Palmer1'2.  Chamber  pressure  varies  in  the  range  2.5-4  bars  while  a  temperature  of  about  1500  K  was 


estimated.  Comparison  of  the  experimental  results  with  the  theoretical  ones  shows:  1)  the  theoretical  mass  flow  rates 
are  lower  than  the  experimental  ones;  2)  the  theoretical  and  experimental  thrust  values  are  in  the  same  range.  This 
point  will  be  discussed  at  the  end  of  this  section. 

In  this  temperature  range  the  kinetics  of  ammonia  presents  large  uncertainties6"9,  even  though  it  has  been  widely 
studied.  The  main  problem  is  that  the  slow  kinetics  prevents  us  from  estimating  the  kinetic  constants  precisely. 
Experimental  works1'2  on  ammonia  thrusters  assume  a  complete  decomposition  of  ammonia  in  the  reservoir,  but  the 
kinetic  scheme  proposed  does  not  lead  to  equilibrium  conditions  in  an  acceptable  time  for  T<2000  K.  Under  these 
conditions  the  most  important  problem  is  that  it  is  possible  that  in  the  inlet  conditions  there  is  not  chemical 
equilibrium  as  usually  happens  in  nozzle  reservoirs.  Because  of  this,  we  also  try  as  an  initial  condition  a  pure 
ammonia  gas. 

Dissociation  inside  the  reservoir  and  the  nozzle  is  very  poor,  so  that  again  in  this  case  the  flow  is  frozen.  As  a 
consequence,  the  flow  properties  depend  only  on  the  Cp/Cv  ratio  and  on  the  mean  molar  mass,  so  the  exit  properties 
are  completely  different  from  the  equilibrium  case.  As  an  example  in  figure  6  we  compare  the  McC  s  computed 
starting  from  an  equilibrium  composition  and  from  pure  ammonia  (non-equilibrium).  It  is  not  possible  to  really 
know  what  composition  prevails  inside  the  reservoir  for  T<2000K,  moreover  some  impurities  or  wall  catalytic 
effects  can  accelerate  the  kinetics,  bringing  the  system  to  equilibrium  faster  than  predicted  by  our  kinetic  model. 

We  try  the  same  comparison  at  T0=2500  K.  In  this  case,  the  mixture  reaches  the  equilibrium  composition  in  the 
reservoir.  Nevertheless,  the  Ma  profiles  in  the  equilibrium  and  non-equilibrium  case  are  again  different  (see  fig.  7). 
This  is  the  demonstration  that  the  history  of  the  gas  in  the  reservoir  is  important. 

Effects  of  initial  conditions  propagate  also  on  the  temperature  profiles.  In  figure  8  we  show  the  normalized 
temperature  profile  calculated  starting  from  an  equilibrium  condition  and  from  an  undissociated  ammonia  (non¬ 
equilibrium)  for  a  reservoir  temperature  of  1500  K  and  a  pressure  of  4  bars.  We  observe  that  the  exit  temperature  is 
much  higher  in  the  case  of  the  equilibrium  conditions.  This  behavior  could  be  due  to  the  lower  mean  molar  mass 
associated  with  the  mixture  in  equilibrium  composition  inside  the  reservoir. 

Different  behaviors  occur  when  the  reservoir  temperature  reaches  2500  K  (see  fig.  9).  In  fact,  at  this  temperature, 
the  reactions  in  the  inlet  are  faster  than  in  the  flow,  therefore  one  can  observe  an  increase  of  the  temperature  due  to 
the  mixture  dissociation.  Even  if  the  ammonia  energy  is  lower  than  that  of  nitrogen  and  hydrogen  molecules,  the 


compression  due  to  the  double  particle  number  density  produces  an  increase  of  the  temperature.  At  the  nozzle  exit, 
the  differences  are  lower  than  in  the  case  of  an  inlet  temperature  of  1500  K. 

For  T0=2500  K  vibrational  kinetics  can  have  appreciable  effects,  therefore  we  perform  our  calculation  including 
the  vibrational  levels.  We  find  similar  behaviors  for  the  relaxation  of  the  vibrational  temperature,  calculated  as  the 
temperature  of  the  first  excited  level  (figure  10).  Starting  from  non-equilibrium  conditions  we  can  observe  an 
increase  of  the  vibrational  temperatures  following  the  gas  temperature.  At  the  nozzle  exit  the  vibrational  temperature 
freezes  at  a  temperature  higher  than  2300  K  just  after  the  throat.  N2  vibrational  temperature  is  practically  frozen. 

As  a  consequence  of  the  previous  results,  we  can  state  that  there  are  some  uncertainties  about  the  inlet 
composition.  To  demonstrate  the  sensitivity  we  have  reported  in  table  7  the  calculated  mass  flow  rate  and  thrust  for 
different  dissociation  degrees  “an”  of  ammonia  for  the  process  NH3  =  1.5  H2  +  0.5  N2  defined  as 

[nh3] 

^  [NH3]+.5[N2] 

that  ranges  from  0  (pure  ammonia)  to  1  (N2  and  Fl2  mixture).  We  note  that  both  quantities  decrease  when  the 
dissociation  degree  increases,  but  the  thrust  changes  more  slowly  than  the  mass  flow  rate,  so  the  specific  impulse 
will  be  higher  for  the  dissociated  case. 

The  last  aspect  that  needs  to  be  investigated  is  the  role  of  internal  degrees  of  freedom,  in  particular  molecular 
vibration,  to  the  measured  quantities  (thrust  and  mass  flow  rate).  The  main  problem  is  that  nothing  is  known  about 
vibrational  relaxation  of  ammonia  molecules.  Therefore  we  compare  the  results  obtained  with  a  frozen  vibrational 
kinetics  with  those  obtained  considering  the  vibrational  states  in  equilibrium  with  the  translational  degrees  of 
freedom  (see  fig.  1 1  and  Tab.  8).  It  can  be  observed  that  the  mass  flow  rate  and  the  specific  impulse  are  weakly 
affected  by  the  internal  state  kinetics  even  if  some  macroscopic  quantities  such  as  Ma  are  strongly  influenced  by  the 
vibrational  contribution,  especially  when  the  inlet  ammonia  concentration  is  high.  This  is  due  to  the  large  number  of 
vibrational  modes  of  ammonia  molecules. 

It  should  be  noted  that  the  new  theoretical  results  with  undissociated  ammonia  (a<r~0)  are  in  better  agreement 
with  the  experimental  mass  flow  rates  even  though  the  comparison  with  thrust  values  is  a  little  worse.  More  detailed 
experimental  data  are  in  any  case  necessary  to  validate  the  present  model,  which  however  at  this  stage  of 
development  represents  a  useful  tool  for  predicting  the  performances  of  ammonia  solar  thrusters.  It  is  also  noted  that 
ref.  1  points  out  evidence  of  receiver  leakage  around  the  seals  during  the  measurements,  which  may  account  for  the 


surprisingly  low  measured  thrust  values.  However,  the  simulation  results  presented  here  provide  a  strong  indication 
that  the  ammonia  propellant  is  in  fact  far  from  completely  dissociated  under  these  conditions. 

D.  High  reservoir  temperature  (T>2500  K) 

In  the  previous  paragraph  we  have  studied  the  behavior  of  the  ammonia  solar  thruster  in  the  working  conditions 
of  the  thruster  investigated  experimentally  in  refs.  1,2.  In  this  paragraph  we  want  to  understand  the  behavior  of  the 
expansion  for  higher  temperature  ranges  (3000  K  <  T0  <5  000  K).  For  temperature  T0>2500  the  ammonia 
dissociation  is  very  fast. 

Moreover  in  the  considered  temperature  range  the  molar  fraction  of  atomic  hydrogen  is  not  negligible.  For  this 
condition,  due  to  the  recombination  of  hydrogen  atoms  and  to  atom-molecule  vibration-translation  collisions,  non¬ 
equilibrium  vibrational  distributions  can  be  observed  and  non-Arrhenius  chemical  rates  are  obtained13,  inducing 
anomalous  behavior  in  the  kinetics. 

Here  we  want  to  focus  on  the  differences  between  the  results  obtained  with  two  different  models: 

a  -  macroscopic  kinetics 

b  -  state  to  state  kinetics. 

In  fig.  12  we  have  compared  the  molar  fractions  of  some  species  calculated  with  the  two  models.  It  is  possible  to 
observe  that  the  H2  molar  fraction  increases  when  the  state  to  state  model  is  considered,  while  the  macroscopic 
model  predicts  a  frozen  concentration.  In  contrast,  NH2  and  NH  molar  fractions  show  smaller  variations  when  the 
state  to  state  approach  is  used. 

This  behavior  affects  the  macroscopic  quantities  such  as  the  gas  speed,  the  temperature  and  the  thrust.  In  fig.  13 
we  have  compared  the  axial  speed  profile  in  the  two  kinetic  models  at  two  different  reservoir  temperatures.  There 
are  appreciable  differences  at  the  nozzle  exit  especially  for  reservoir  temperature  T0=5000  K.  This  effect  is  the 
consequence  of  the  stronger  recombination  of  hydrogen  atoms  producing  H2  molecules  as  can  be  observed  in  figure 
12.  Similar  behaviors  can  be  observed  in  the  gas  temperature  profile  (fig.  14)  and  in  the  thrust  (fig.  15).  An 
important  result  is  that  the  calculated  thrust  (and  therefore  the  performance)  is  higher  when  the  state  to  state 
vibrational  kinetic  model  is  included  in  the  calculation. 

As  for  the  low  temperature  case,  the  molar  fraction  of  NHX  species  is  small  so  that  their  influence  to  the  flow 
properties  is  negligible.  However,  non-equilibrium  vibrational  distributions  of  such  molecules  could  affect  the  gas 
composition  due  to  deviations  of  the  global  rates  from  the  Arrhenius  trend. 


IV.  Conclusion 


In  this  paper  we  have  investigated  the  chemical  kinetic  aspects  of  an  ammonia  mixture  for  application  to  solar 
thermal  thrusters.  We  have  obtained  results  that  contradict  the  common  assumption  of  considering  inlet  conditions 
in  equilibrium.  If  the  kinetics  of  the  mixture  are  very  slow,  as  in  the  case  of  ammonia,  the  gas  composition  can  be  far 
from  equilibrium.  The  comparison  of  measured  mass  flow  rates  with  a  range  of  assumed  dissociation  fractions  also 
supports  the  conclusion  that  ammonia  propellant  at  this  temperature  may  be  substantially  undissociated.  At  low 
reservoir  temperature  (T0<2500  K),  that  is  the  working  condition  for  the  thruster  in  the  laboratory,  the  kinetics  of 
internal  states  is  negligible,  and  the  macroscopic  and  state  to  state  models  gives  the  same  results,  mainly  because  the 
flow  is  almost  frozen.  On  the  other  end,  at  higher  reservoir  temperature  (T0>3000  K),  when  the  hydrogen 
dissociation  is  important,  the  vibrational  kinetics  becomes  very  important.  In  these  conditions  the  macroscopic  and 
the  state  to  state  model  predict  different  results,  mainly  due  to  the  different  atomic  recombination. 

We  must  emphasize  that  this  study  is  preliminary.  Some  properties  of  the  mixtures  have  been  neglected,  such  as 
the  vibrational  kinetics  of  the  NHX  molecules  and  the  formation  of  NH  molecules  by  vibrationally  excited  N2  and  H2. 
The  model  can  be  improved  including  a  multitemperature  approach  to  account  for  the  vibration  of  polyatomic 
molecules  and  by  adding  new  state  to  state  reactions  between  N2  and  H2. 
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Table  1  Species  inserted  and  their  physical  properties4  n.  Rot  is  the  number  of  independent 
rotational  axis,  Vib  the  number  of  vibrational  modes  and  Ef  the  formation  energy. 


Species 

Mass  (AMU) 

Rot 

Vib 

Ef(eV) 

nh3 

17 

3 

6 

-0.403242 

n2h3 

31 

3 

9 

4.30994 

nh2 

16 

3 

3 

2.00285 

NH 

15 

2 

1 

3.90225 

n2 

28 

2 

1 

0.00000 

h2 

2 

2 

1 

0.00000 

N 

14 

0 

0 

4.87950 

H 

1 

0 

0 

2.23910 

Table  2  Spectroscopic  constants  (in  cm'1)  used  in  the 
calculation  of  the  vibrational  level  energies5  (see  eq.  1). 


Parameters 

n2 

H  2 

®e 

2358.57 

4464.7 

roexe 

14.324 

146.93 

caeye 

-0.00226 

4.7108 

COeZe 

0.0 

-0.22712 

Ediss 

9.7639 

4.4772 

V  max 

48 

15 

Table  3  Arrhenius  coefficients  for  the  reactions  inserted  in  the  model6'11.  K0  has  the  dimension  of 
cm3(n  l)/s,  where  n  is  the  number  of  particles  involved  in  the  reaction,  Ea  is  expressed  in  K  and  8  is 

a  dimensionless  parameter  (eq.  2). 


process 

K0 

5 

Ea 

Ref. 

pi 

1.53e-8 

0 

42400 

6 

p2 

1.00e-9 

0 

47000 

4 

p3 

3.20e-13 

0.67 

1720 

5 

p4 

1.00e-32 

0 

0 

7 

p5 

1.00e-34 

0 

0 

7 

p6 

2.30e-13 

0.67 

2160 

5 

P7 

6.64e-12 

0 

2800 

6 

p8 

1.00e-33 

0 

0 

7 

p9 

1.70e-12 

0.68 

957 

5 

plO 

1.77e-ll 

0 

0 

5 

pll 

6.60e-13 

0.55 

957 

5 

pl2 

1.34e-31 

-0.60 

0 

5 

pl3 

2.68e-31 

-0.60 

0 

5 

pl4 

7.44e-32 

-0.50 

0 

5 

pl5 

3.31e-27 

-1.50 

0 

5 

pl6 

1.33e-12 

0.5 

10850 

4 

P 1 7 

1.70e-8 

0. 

21000 

4 

p  1 8 

2.60e-12 

0 

0 

4 

p!9 

1.70e-13 

0 

0 

4 

Table  4  Coefficients  used  to  calculate  the  formation  equilibrium  constants  when  the  pressures 


are  expressed  in  Pa  (eq.  4). 


Species 

K,, 

KP 

Op 

Ke 

9e 

nh3 

51.8293 

-13.1631 

0.960208 

-7.97367 

628.095 

nh2 

19.6666 

63.4542 

1.01863 

13.5452 

311.287 

NH 

-5.14978 

85.2911 

0.998065 

-2.21874 

423.673 

N 

28.3772 

-114.298 

0.997620 

-2.14877 

4371.66 

H 

26.1506 

-52.7759 

0.996362 

-2.71911 

2157.34 

Table  5  Computed  mass  flow  rate  and  thrust  of  the  nozzle  as  a 
function  of  the  inlet  pressure.  The  inlet  temperature  is  T0=1500  K 
except  when  marked  with  the  “*”  symbol  assigned  to  T0=2500  K.  At 
the  inlet  we  suppose  an  equilibrium  composition. 


Po(bars) 

mass  flow  rate  (g/s) 

thrust  (mN) 

1 

2.1374e-2 

63.379 

2 

4.2766e-2 

126.71 

3 

6.4122e-2 

189.99 

4 

8.5496e-2 

253.58 

4* 

6.6917e-2 

256.05 

Table.  6  Mass  flow  rate  and  thrust  of  the  nozzle 


determined  experimentally  in  different  test  cases. 


mass  flow  rate  (g/s) 

thrust  (mN) 

.140 

219 

.147 

223 

.117 

183 

.121 

190 

Table  7  Computed  mass  flow  rate  and  thrust  as  a  function  of  the  inlet  dissociation 
degree  (aD)  of  the  ammonia.  Inlet  pressure  and  temperature  are  P0=4  bars  and 
T0=1500  K.  Vibrational  degrees  of  freedom  are  frozen. 


aD 

mass  flow  rate  (g/s) 

thrust  (mN) 

ISp 

0.0 

1.2120e-l 

261.56 

220 

0.2 

1.1099e-l 

260.72 

239 

0.4 

1.0317e-l 

259.64 

256 

0.8 

9.1461e-2 

257.96 

287 

1.0 

8.5496e-2 

253.52 

302 

Table  8  Computed  mass  flow  rate  and  thrust  as  a  function  of  the  inlet  dissociation 


degree  (aD)  of  the  ammonia.  Inlet  pressure  and  temperature  are  P0=4  bars  and 
To=1500  K.  Vibrational  degrees  of  freedom  are  in  equilibrium  with  the 

translational  ones. 


aD 

mass  flow  rate  (g/s) 

thmst  (mN) 

ISp 

0.0 

1 . 1 364e- 1 

282.55 

253 

0.2 

1.0438e-l 

280.42 

274 

0.4 

9.7589e-2 

276.91 

289 

0.8 

8.7962e-2 

269.78 

312 

1.0 


8.2806e-2 


261.59 


322 


radius  (m) 


inlet  throat 


x(m) 


Figure  1.  Nozzle  profile. 
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Figure  5:  Pressure  profile  along  the  nozzle  as  a  function  of  the  inlet  pressure  for  T0=lSOO  K  and  starting 


from  an  inlet  equilibrium  composition. 
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Figure  6:  Ma  profile  compared  for  inlet  equilibrium  and  pure  ammonia  (non-equilibrium)  case  (pressure  and 


temperature  conditions  as  in  fig  2). 
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Figure  7:  Same  as  in  fig  6  with  reservoir  temperature  T0=2SOO  K  and  pressure  Po=4  bars. 
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Figure  8:  Reduced  temperature  profile  compared  for  inlet  equilibrium  or  pure  ammonia  (non-equilibrium) 
case  (pressure  and  temperature  conditions  as  in  fig  2). 
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Figure  9:  Reduced  temperature  profile  compared  for  inlet  equilibrium  or  pure  ammonia  (non-equilibrium) 
case,  for  a  reservoir  temperature  T0  =  2500K  and  pressure  Po=4  bars.. 
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Figure  11:  Ma  profile  compared  for  inlet  equilibrium  and  pure  ammonia  (non-equilibrium)  case  (pressure 
and  temperature  conditions  as  in  fig  2).  The  vibrational  degrees  of  freedom  are  in  equilibrium  with  the 


translational  ones. 
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Figure  13:  Velocity  profiles  calculated  with  different  kinetic  models  and  inlet  temperature. 
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Figure  14:  Temperature  profiles  calculated  with  diffe 
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Figure  15:  Thrust  profiles  calculated  with  different  kinetic  models  and  inlet  temperature. 


